We extend the core-accretion model of giant gaseous planets by Pollack et al. f fToM] ! 
to include migration, disc evolution and gap formation. Starting with a core of a fraction 
of an Earth's mass located at 8 AU, we end our simulation with the onset of runaway 
gas accretion when the planet is at 5.5 AU 1 Myr later. This timescale is about a factor 
ten shorter than the one found by Pollack et al. (199^) even though the disc was less 
massive initially and viscously evolving. Other initial conditions can lead to even shorter 
timescales. The reason for this speed-up is found to result from the fact that a moving 
planet does not deplete its feeding zone to the extend of a static planet. Thus, the 
uncomfortably long formation timescale associated with the core-accretion scenario can 
be considerably reduced and brought in much better agreement with the typical disc 
lifetimes inferred from observations of young circumstellar discs. 
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1 Introduction 

The current paradigm for the formation of giant gaseous planets is based on the so-called 
core accretion model in which a growing solid core reaches a critical mass and accretes 
rapidly a massive atmosphere (Pollack et al. jl996j . hereafter referred to as P96). While 
this model has many appealing features, it suffers at least from two shortcomings which, 
as we shall show later, are actually coupled. 

First, the timescale (close to 10 Myr) found by P96 to form Jupiter at its present 
location is uncomfortably close to the typical lifetime of protoplanetary discs which is 
believed to be of the order of 1-10 Myr (Haisch et al. 12001^ . This timescale problem 
has led others to look for more rapid formation mechanisms based on direct gravitational 
collapse (Boss j2002j ). Second, P96 assumed that the giant planets of our solar system have 
been formed where they are observed today. However, the discovery over the last decade 
of extrasolar planets at very short distances to their parent star, has open the possibility 
that planets may actually migrate over large distances (Lin et al. jl99()j : Trilling et al. 
jl998j : Papaloizou & Terquem jl999j . hereafter referred to as PT99). The timescale of 
migration is still very uncertain, but conservative estimates give values between 0.1 and 
10 Myr. This timescale is therefore comparable to the planet formation timescale and the 
disc lifetime and thus migration cannot be neglected in a self-consistent picture. 

In this Letter, we extend the core accretion model of P96 to include migration, disc 
evolution and gap formation and show that the formation timescales of giant planets can 
be reduced by factors of ten or more. 

2 The Model 

Our model consists of three different modules that calculate: 1) the disc structure and 
time evolution, 2) the interaction of planetesimals with the atmosphere of the planet, and 
3) the internal structure of the planet. A detailed description of the model and of the 
tests performed will be given in a forthcoming paper, here we only briefly point out the 
main features of the different modules. 
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2.1 Disk structure and evolution 

This module determines the structure (both vertical and radial) of a protoplanetary disc 
modelled as an a— disc. For this, we follow the method described in PT99, and refer the 
reader to this paper for more details. Using their procedure we calculate, as a function 
of the surface density S and the radial distance to the star r, the mid-plane temperature 
and pressure Tmid{r,T,), the mean viscosity i>(r, S) and the disc density scale 

height H{r,T,). 

The time evolution of the disc is governed by a diffusion equation, modified to take into 
account the momentum transfer between the planet and the disc. The rate of momentum 
transfer between the planet and the disc is calculated using the formula derived by Lin &: 
Papaloizou ( |1986j ) (equation 14). 

2.2 Migration 

Gravitational interactions between the growing protoplanet and the disc lead to inward 
migration and possibly gap formation (Lin & Papaloizou |1986j . Ward |1997j . Tanaka et 
al. ^2002,). For low mass planets, the migration rate is linear with mass (type I migration. 
Ward |1997j ). Higher mass planets open a gap and the migration rate is set by the viscosity 
independently of planetary mass (type II migration. Ward I1997j ). 

While the general physical understanding of the origin of migration is clear, the actual 
migration rates obtained for type I migration especially are so short that all planets should 
actually be destroyed by the central star long before the disappearance of the gaseous disc. 
Tanaka et al. ( j2()()2j ) have performed new analytical calculations of type I migration, in 
two or three dimensional discs and found longer migration timescales but unfortunately 
still too short to ensure survival. Indications for longer type I migration timescales can 
be found in calculations by Nelson & Papaloizou ([2003 ). As suggested by these authors, 
torques exerted on at least small mass planets (Mpianet < 3OM0) embedded in turbulent 
MHD discs are strongly fluctuating resulting in a slow down of the net inward motion. 

These considerations seem to indicate that the actual migration timescales may in 
fact be considerably longer than originally estimated by Ward ( I1997j ) or even by Tanaka 
et al. ( |2002j ). For these reasons, and for lack of better knowledge, we allow us the 
freedom to use for type I migration the formula derived by Tanaka et al. ( j20()2j ) reduced 
by an arbitrary numerical factor // (set to 1/30 in this paper). Tests have shown that 
provided this factor is small enough to allow planet survival, its actual value does not 
change our main conclusion but only the extend of migration. For type II migration, the 
inward velocity is set by the viscosity of the disc. Migration type switches when the planet 
becomes massive enough to open a gap in the disc. In our work, this transition occurs 
when the Hill radius of the planet becomes greater than the density scale height H of the 
disc. 
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2.3 The planetesimals 

2.3.1 Interaction with the growing atmosphere 

In our second module we determine the trajectory, the energy and mass loss of planetes- 
imals falling through the atmosphere of the planet under the influence of gravity and 
aerodynamic drag forces. The drag coefficient is calculated (assuming a sphere) as a func- 
tion of the local Mach and Reynolds number using the equations given by Henderson 
f |1976j ^. The loss in kinetic energy results in a local heating of the planet's atmosphere 
which enters in the calculation of the internal structure. Given the size of the planetesimals 
considered here (100 km), ablation is found to be negligible and deposition of mass occurs 
almost entirely due to fragmentation which occurs when the pressure at the stagnation 
point becomes larger than its tensile strength. For the envelopes considered here, this 
occurs typically at or quite close to the core making the difference between the sinking 
and the no-sinking case as defined by P96 much smaller. Note that in the current version 
of the code we do not include the effect of heavy elements enrichment in the opacity and 
the equation of state. 

2.3.2 Accretion rate of soUds 

We assume that due to the scattering effect of the planet, the surface density of planetes- 
imals is constant within the current feeding zone but decreases with time proportionally 
to the mass accreted by the planet. The feeding zone is assumed to extend to a distance 
Qf = 4 X i?H on each side of the planetary orbit, where -Rh = (-^planet ) Opianet 
is the Hill radius of the planet, Cpianet the star-planet distance and the mass of the 
central star. We use the expressions in the Appendix B and C of Greenzweig & Lissauer 
f |1992j ^ to calculate the gravitational enhancement factor Fg. The solid accretion rate is 
given by: = R'^J^pflFg/ (27r) where Sp is the surface density of solids at r = Opianet 

and Rc is the capture radius of the planet, calculated as in P96, using our second module. 
For the inclination and eccentricities of the planetesimals we adopt the same values as in 
P96. 

2.4 Protoplanet structure and evolution 
2.4.1 Internal structure equations 

In the third module, we calculate the internal structure of the planet including a growing 
core (and the inner luminosity due to the accretion of planetesimals) and gas accretion due 
to both the contraction of the envelope and the increase of the outer radius of the planet 
(see P96). The standard equations of planet evolution are solved, using opacities from 
Bell & Lin (1994) and the equation of state (EOS) from Chabrier et al. (1992). Note that 
in the present models we neglect the change of entropy while computing the luminosity 
gradient. This approximation is justified by the fact that until runaway gas accretion, the 
energy deposition by infalling planetesimals largely dominates the energy budget. 

The core radius is set to -Rcore = (3Mcore/47rpcore)^'^^ and the core luminosity is equal 
to the remaining energy of planetesimals after having passed through the atmosphere. The 
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density of the core is fixed to 3.2g/cm^. The two outer boundary conditions are given by 
requiring that the disc and the planet joint smoothly at the outer radius, i.e. Pgurf = Pmid, 
and Tsurf = Tmid where Pmid and Tmid are the local mid-plane values obtained from the 
disc structure calculations (see above). These two conditions are valid as long as the disc 
can supply enough mass to keep the outer radius equal the Hill (or the accretion) radius 
(see next section). 

2.4.2 Gas accretion rate 

The gas accretion rate onto the planet is determined by the condition: i?pianet = min(i?H, -Racer) 
where i?pianet is the outer radius of the planet, and -Racer = CMpianet /C's j Cs is the local 
sound velocity in the disc (see P96). At each timestep, we iterate on the mass of the 
envelope (and then the total mass) until this condition is met. 

In reality, the condition -Rpianet = niin(-RH; -Raeer) can only be satisfied if the disc can 
actually supply enough gas to the planet. When a gap opens, the gas available for the 
planet to accrete drops significantly. Hence, we limit the gas delivery rate from the disc 
to the planet to the maximum value given by Veras & Armitage ( j2UU3j ). 

2.5 Tests of the model 

Various tests have been performed to validate our entire model and we will describe them 
in detail in a forthcoming paper. Let us here only mention that our disc structure calcu- 
lations reproduce the results given in PT99 whereas the module dealing with the internal 
structure of the planet is in agreement with the results of Bodenheimer & Pollack ( jl986j ). 
Our planetesimal accretion module has been tested extensively by comparing results of 
simulations of impacts on Earth (Hills & Goda jl99,Sj ). on Venus (Zahnle jl992j ). and 
Jupiter (eg. comet Shoemaker-Levy 9, MacLow &; Zahnle jl994j '). Moreover, the entire 
code has been tested using the same initial conditions as P96 (case Jl) turning disc evo- 
lution and migration off. In this case, we obtain a formation time of 9 Myr, in close 
agreement to their result. 

3 Results 

3.1 Initial conditions 

We start with a solid core of 0.6M^ located at 8 AU which will eventually lead after ~ 1 
Myr to a giant planet now located around 5.5 AU but still migrating inward. 

To allow an easy comparison with P96, we chose an initial disc surface density profile 
given by a power law S oc r~^, where r is the distance to the star, the constant being 
chosen to yield S = 525g/cm^ at 5.2 AU. The viscosity parameter a is set to 2 x 10~^ 
which yields a typical evolution time of the disc of 4 Myr. The gas-to-dust ratio is equal 
to 70. Our initial disk therefore corresponds closely to case J2 of P96 for which they found 
a formation time of about 50 Myr. The numerical parameters are // (reduction of type I 
migration) equal to 1/30 and /a, the numerical factor in the expression of the momentum 
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Figure 1: Total mass of heavy elements (core+envelope) and mass of the envelope (H/He) 
as a function of time, for three models. Solid lines: migrating case, without gap formation 
(case 1); dashed lines: migrating case, with gap formation (case 2); and dot-dashed lines: 
non migrating case, without disk evolution nor gap formation (case 3). The heavy line 
gives the distance to the central star as a function of time (case 1). The kink at 0.8 Myr 
is the transition from type I to type II migration. 

transfer between the planet and the disc, set to 0.05. This latter choice gives a transition 
from type I to type II migration when the reduction of S due to momentum transfer is 
around 30%. We have checked that our conclusions remain valid if /a is set to (no gap 
formation) . 

3.2 Formation timescales 

Figure n shows the mass of heavy elements and the mass of H/He accreted by the planet 
as a function of time. Note that the mass of heavy elements does not correspond to the 
core mass since some fraction of the planetesimals are being destroyed while traversing 
the envelope and never reach the core. 

As in P96, the formation timescale is essentially determined by the time necessary to 
reach the runaway accretion phase which occurs shortly after the crossover mass (mass of 
core equals mass of envelope), Mcross; has been reached. In the case where migration and 
disk evolution have been switched off, this occurs after ~ 31 Myr, compared to 48 Myr 
in P96. This difference provides a measure for the sensitivity of the results to differences 
in physical and numerical approximations used in both approaches. In particular, we use 
different EOS and, more importantly, opacity laws (which has been shown to have a huge 
influence on formation timescale, see model J6 of P96). To properly derive the effect of 
migration and disk evolution, we will compare the formation timescales obtained with our 
own value (~ 30 Myr). 

Allowing for migration and disc evolution, we obtain a formation time of about ~ 1 
Myr, i.e. thirty times faster than in our identical model in which migration and disc 
evolution has been switched off. This is also an order of magnitude faster than the pre- 
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Figure 2: Solid and gas surface densities for two of our simulations. Solid line: initial gas 
surface density; dotted line: solid surface density for case 3 after 1.2 Myr; dashed line: gas 
surface density at the end of case 2; dot-dashed line: solid surface density for the same 
model, at the same time. The solid surface densities are multiplied by 70. The dot gives 
the position of the planet. 

ferred model (Jl, S normalized to 700g/cm^) of P96 even though their initial disc was 
significantly more massive than ours and not viscously evolving. The main reason for this 
speed-up is that owing to migration the feeding zone is not as severely depleted as in P96. 
In their models, it is this depletion which was responsible for the long time needed by the 
core to reach critical mass and start runaway gas accretion. In our models, the moving 
planet always encounters new planetesimals and thus its feeding zone is never emptied. 
To illustrate this important point, we show in Figure 121 the initial and final disc profiles 
(for both the gaseous and the solid component). 

We note that the crossover mass obtained in our calculation (-^ 22Mq)) is significantly 
larger ^ than the one obtained by P96 (case Jl and J2, respectively IGM^ and ~ 
IO.5M0). This is to be expected since the crossover mass is a growing function of the 
solid accretion rate. Migration, besides preventing the depletion of the feeding zone, also 
causes higher accretion rates (and then higher crossover masses) than the corresponding 
no-migrating case. However, we note that this crossover mass is significantly lower than 
the one obtained by P96 in their model J3 (S normalized to 1050g/cm^, Mcross ~ 29M^), 
which started with a disk twice as massive as ours in order to obtain formation timescales 
comparable to ours. 

4 Summary and discussion 

We have calculated the formation of giant planets up to runaway gas accretion and studied 
the effect of migration and gap formation on the resulting formation timescale. Our main 

^However, the final mass of heavy elements is still compatible with the one derived by Guillot f |1999| 'l. 
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result is that phase 2 described in P96 is suppressed, and this leads to a formation time 
around 1 Myr for the initial conditions chosen here. This is roughly an order of magnitude 
faster that the favored model of P96 even though our initial disc had less mass than the 
one used by P96. 

We have performed many tests and convinced ourselves that this speed-up is robust 
against changes in various parameters. For example, in a calculation in which the reduction 
of type I migration (//) is set to 0.1, an embryo starting at 15AU undergoes runaway 
accretion in less than 3 Myr. The assumed size of the planetesimals plays a critical role 
as already noted by P96. For example, assuming planetesimals of 10km instead of 100km 
leads to runaway accretion after only 0.3 Myr! As already mentioned in section !^. 3. 11 since 
the effect of ablation is negligible in our calculations, mass loss of planetesimals occurs 
very deep in the envelope. The switch from sinking to no-sinking approximation (see 
P96) has then a small effect. In the migrating case (without gap formation) we obtain a 
formation timescale of ~ 1.2 Myr using the sinking approximation, compared to ~ 1 Myr 
in the no-sinking case. 

Thus, migration, besides explaining the presence of giant planets at short distances 
to their stars, also plays an important role in the formation process itself. By ensuring 
that the feeding zone is never depleted, migration suppresses almost entirely the pro- 
tracted phase 2 therefore shortening enormously the formation time. This of course, does 
not preclude other effects such as reduction of opacity or formation of vortices prior to 
planet formation to further reduce this timescale. The formation of giant planets through 
the core-accretion scenario can therefore proceed over timescales in good agreement with 
disk lifetimes, and without having to consider disks significantly more massive than the 
minimum mass solar nebula. 

Finally, we point out that the formation of giant planets appears only possible if the 
currently available type I migration rates are reduced by at least a factor of 10. This 
suggests that their might still be a serious problem in our understanding of this type of 
migration. 
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